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' Abstract. We report on the first near-IR [Fe II] line imaging survey of extragalactic SNRs. Observations of a sample 

^ of 42 objects drawn from an optically-selected catalogue of SNRs in M33 provide evidence for a wide range in the 

QQ , [Fe II] A1.644 fim luminosities. This can be understood as being primarily due to variations in the chemical abundances 

■ and density of the local ISM, although shock conditions may also play a significant role. We briefly discuss how these 
' results may be used to better calibrate the supernova rate of star-forming galaxies. 

o _ 

^ . 1 Why this survey? 

The use of the strong near-IR [Fe II] lines as a mean of estimating the supernova rate of starburst 
^ ■ galaxies has become increasingly popular over the last decade (e.g., Calzetti 1997). However, one 

^ ', serious limitation of this technique lies in the fact that one has to adopt typical values for the [Fe II] 

' luminosity and [Fe II]-emitting lifetime of a single supernova remnant (SNR), whereas a wide span 

of values are observed. For instance, the [Fe II] luminosities of SNRs in the starburst galaxy M82 
(Greenhouse et al. 1997) are two orders of magnitude higher than in "quiescent" galaxies, such as 
our Galaxy or the LJVIC (Oliva et al. 1989). One order of magnitude higher values are also observed 

■ in NGC 253 (Forbes et al. 1993). Therefore, before applying this technique to star-forming galaxies 
that presumably exhibit widely different properties (e.g., metallicity) , one must first understand what 
parameters control the [Fe II] properties of individual SNRs. In an effort to address this issue, we 
report here on preliminary results of a [Fe II] A1.644 /im line imaging survey of an heterogeneous 
sample of SNRs in M33. 

2 Observations 

Our observations were obtained during two observing runs at the Canada-France-Hawaii telescope 
(CFHT) in 1997 and 1998. Two different infrared cameras were used: MONICA in 1997 and REDEYE 
in 1998. A total of 42 objects drawn from the optically-selected atlas of Gordon et al. (1998) has been 
observed. Particular attention has been paid to observe SNRs with widely different properties. The 
standard reduction procedure for infrared imaging has been applied (see, e.g., Hodapp et al. 1992). 

3 Relationship between the [Fe II] emission and other SNR properties 

The broad range of [Fe II] A1.644 fim luminosities is made clear by the fact that — while several SNRs 
present an outstanding level of [Fe II] emission (see Fig.l) — many of them remain undetected in our 
survey. 
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Fig. 1. Comparison of the H {left panel) and [Fe II] A1.644 fim (right panel) images of SNR 28 (identification from 
Gordon et al. 1998). The SNR is indicated by an arrow. The field of view is 20" X 20". North is up, and east is to the 
left. The integration times are 300 s and 1 500 s for H and [Fe II] , respectively. 



The [Fe II] A1.644 /im luminosities (or upper limits) are presented in Table 1. It can be seen that a 
(nearly) two orders of magnitude range is observed in these values. Lumsden & Puxley (1995) reported 
near-IR spectroscopic observations of 9 of our program targets. There is a good agreement between 
the two estimates.0 

In order to pin down the physical processes that produce the [Fe II] emission, we sought for 
correlations between the [Fe II] A1.644 /im luminosities and various quantities intrinsic to the SNRs, 
such as their age or optical-line properties, for instance. We made use of the generalized Kendall's tau 
correlation technique which allows us to investigate statistical correlations between datasets including 
upper limits (Isobe, Feigelson, & Nelson 1986). Table 2 summarizes the results of this statistical 
analysis for all quantities under consideration, along with the source of these data. 

In a preliminary effort to calculate the dynamical age of the SNRs, we have simply assumed that 
they all undergo Sedov- Taylor expansion, as suggested by Gordon et al. (1998). Without a detailed 
knowledge of the evolutionary status, initial blast energy of the SNR, and of the ambient gas density 
in which it evolves, age estimates are notoriously uncertain. Therefore, the complete lack of any 
relationship between the level of [Fe II] emission and the dynamical age of the remnant (see Table 2) 
should be treated at this stage of the data analysis with caution. 

Our data support the idea that environmental rather than evolutionary effects are more 
important in controlling the [Fe II] properties. A tight correlation is found between the electronic 
densities derived from the [S II] AA6717, 6731 doublet and L{[Fc II] A1.644 ^m). This confirms that the 
[Fe II] emission arises in the dense postshock region. The brightest [Fc II] remnants are characterized 
by ne([S II]) « 10"^ cm^'', which translates to ne{[Fe II]) « lO"^-^ cm^"^ (e.g., Lumsden & Puxley 
1995). This requires preshock densities significantly higher than canonical values for the ISIVI (rig « 
0.1-10 cm~'^), approaching values typical of molecular material. Radio emission from SNRs is due to 
synchrotron emission, and is thus primarily a function of the number of available electrons and of 
the magnetic field strength. The tendency for the objects with high radio fluxes to be strong [Fe II] 
emitters may thus be secondary in nature, and merely reflect the fact that dense environments both 
enhance the [Fe II] and radio emission (provided that the magnetic properties are roughly similar 
within the sample). There is indeed some indication of a correlation between and 6*20 (Table 2). 

A tight correlation is also observed between L([Fe II] A1.644 /im) and the generalized "metal 



^Their [Fe II] A1.257 fim fluxes have been converted into [Fc II] A1.644 ^m fluxes assuming a ratio: [Fe II] A1.257 
/im/[Fe II] A1.644 fim = 1.36 (Nussbaumcr & Storey 1988). 
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70.6 ± 28.6 
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< 432 
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< 92.2 
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< 225 
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< 19.5 
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< 616 
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< 1040 


80 


< 306 
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< 230 
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< 244 


87 


< 103 


94 


189 ± 51 
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Tab. 1. [Fe II] A1.644 /im luminosities of SNRs in M33.° 



" Assuming a distance to M33 of 840 kpc (Freedman, Wilson, & Madore 1991). The upper limits are 
quoted for a statistical significance of 3a. The error bars for the detected objects are la uncertainties. 
^ Identifications from Gordon et al. (1998). 



abundance", A, as defined by Dopita et al. (1984). This most likely indicates that [Fe II] emission 
is enhanced in regions characterized by high ISM chemical abundances. Because of the substantial 
metallicity gradient of M33 (Garnctt ct al. 1997), we should therefore observe a strong negative 
correlation between L{[Fe II] A1.644 /jin) and the SNR galactocentric distance. However, such evidence 
is lacking in our data (false alarm probability of 8.5 %), probably owing to the large scatter in the 
GCD-metallicity relation (Smith et al. 1993). Although iron-bearing grains are easily destroyed by 
shocks (Draine & Salpeter 1979), both theory and observations suggest that the efficiency of this 
process may vary from one SNR to another (Jones, Tielens, & HoUenbach 1996; Oliva et al. 1999). 
The observed correlation between A and L{[Fe II] A1.644 /xm) may also suggest that [Fe II] emission 
is favoured by strong shocks that return a large fraction of the dust grains into the gas phase. In this 
case, we should expect a correlation between the [Fe II] A1.644 ^m luminosities and the expansion 
velocity Vbuik (to first order proportional to the shock speed). This correlation is also not clear in our 
data, possibly because our range of observed Vtuik is not broad enough. Different electronic densities 
can also introduce noise in the correlation. 

Although [Fe II] A1.644 /xm is believed to probe somewhat cooler, denser material than the optical 
[S II] doublet (e.g., Oliva et al. 1989), the correlation between the [Fe II] A1.644 ^m and [S II] fluxes 
suggests that the line-formation regions of these two species are not drastically different and that 
similar physical processes give rise to both emissions. All optical-line fluxes correlate extremely well 
with L{[Fe II] A1.644 /xm). 

4 Discussion 

In the light of the results presented above, we conclude that density and (possibly to a less extent) 
abundance effects arc likely to play a pivotal role in controlling the [Fe II] emission in SNRs. This 
may help to understand the differences that exist between the near-IR [Fe II] properties of SNRs in 
galaxies with moderate star- formation activity (such as M33) and star burst galaxies. While the present 
study demonstrates that the brightest remnants have luminosities of the order of 700 Lq, SNRs in the 
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Tab. 2. Results of the statistical analysis. N: Number of points included in the calculations. 'P{X, Y): Probability for a 
chance correlation between the variables X and Y. Note that the objects at large galactocentric distances, GCD, show 

a tendency for lower L{\Fe II] A1.644 ^m) values. 

X Y N V{X, Y) 



A1.644 /zm) 


GCD" 


42 


8.5 % 


II 


He" 


34 


0.1 % 


II 


Dynamical age" 


42 


82 % 


II 


Metal abundance'' 


20 


0.2 % 


II 


Vhulk" 


8 


8.3 % 


II 


T ([S II] A6717)'^ 


34 


0.1% 


II 




35 


1.5 % 


II 


520*= 


35 


2.2 % 




^20^ 


34 


5.6 % 



Prom Gordon et al. (1998). The dynamical ages have been derived using their quoted optical 

diameters. 

^ From Blair & Kirshner (1985) and Smith et al. (1993). 
Prom Blair, Chu, & Kennicutt (1988). 

From Smith ct al. (1993), Gordon et al. (1998), and Blair & Kirshner (1985). 
^: From Duric et al. (1993) and Gordon et al. (1999). 



starburst galaxy M82 have luminosities up to 1.6 x 10^ Lq (Greenhouse ct al. 1991, 1997). Although 
it remains to be seen whether the two samples trace two similar populations of SNRs, we propose that 
this dichotomy may be explained by the dense, dusty environments prevailing in the nuclear regions 
of these galaxies. In particular, it is conceivable that processing (via grain destruction) of the large 
reservoir of dust grains greatly enhances the [Pe II] emission with respect to more dust-free galaxies. 



References 

Blair, W. P., Chu, Y. -H., Kennicutt, R. C. (1988), in Supernova Remnants and the Interstellar Medium, eds. R. S. 

Roger & T. L. Landecker, Proc lAU CoUoq. 101, Cambridge Univ. Press, Cambridge, p. 193 
Blair, W. P., & Kirshner, R. P. (1985) Astrophys. J. 289, 582. 
Calzetti, D. (1997) Astron. J. 113, 162 

Dopita, M. A., Binette, L., D'Odorico, S., & Bcnvcnuti, P. (1984) Astrophys. J. 276, 653 
Draine, B. T., & Salpeter, E. E. (1979) Astrophys. J. 231, 438 

Duric, N., Viallcfond, F., Goss, W. M., & van dcr Hulst, J. M. (1993) Astron. Astrophys., Suppl. Ser. 99, 217 
Forbes, D. A., Ward, M. J., Rotatiuc, V., Blietz, M., Genzel, R., Drapatz, S., van der Werf, P. P, & Krabbe, A. (1993) 

Astrophys. J., Lett. 406, Lll 
Freedman, W. L., Wilson, C. D., & Madore, B. F. (1991) Astrophys. J. 372, 455 

Garnett, D. R., Shields, G. A., Skillman, E. D., Sagan, S. P., & Dufour, R. J. (1997) Astrophys. J. 489, 63 
Gordon, S. M., Kirshner, R. P., Long, K. S., Blair, W. P., Duric, N., & Smith, R. C. (1998) Astrophys. .J., Suppl. Ser. 
117, 89 

Gordon, S. M., Duric, N., Kirshner, R. P., Goss, W. M., & Viallefond, F. (1999) Astrophys. J., Suppl. Ser. 120, 247 
Greenhouse, M. A., Woodward, C. E., Thronson Jr., H. A., Rudy, R. J., Rossano, G. S., Erwin, P., & Puetter, R. C. 

(1991) Astrophys. J. 383, 164 
Greenhouse, M. A., et al. (1997) Astrophys. J. 476, 105 

Hodapp, K.-W., Rayner, J., & Irwin, E. (1992) Publ. Astron. Soc. Pac. 104, 441 

Isobe, T., Feigelson, E. D., & Nelson, P. I. (1986) Astrophys. J. 306, 490 

Jones, A. P., Tielens, A. G. G. M., & HoUenbach, D. J. (1996) Astrophys. J. 469, 740 

Long, K. S., Charles, P. A., Blair, W. P., & Gordon, S. M. (1996) Astrophys. J. 466, 750 

Lumsden, S. L., & Puxley, P. J. (1995) Mon. Not. R. Astron. Soc. 276, 723 

Nussbaumer, H., & Storey, P. J. (1988) Astron. Astrophys. 193, 327 

Oliva, E., Moorwood, A. F. M., & Danziger, I. J. (1989) Astron. Astrophys. 214, 307 

Oliva, E., Lutz, D., Drapatz, S., & Moorwood, A. F. M. (1999) Astron. Astrophys., Lett. 341, L75 

Smith, R. C, Kirshner, R. P., Blair, W. P., Long, K. S., & Winkler, P. F. (1993) Astrophys. J. 407, 564 



